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Abstract 

We present measurements of time-dependent CP asymmetries in (/)(1020)i^*^, rj'K^, 

K^K^K^, A'OvrO, /o(980)A'^, uj{782)K'^ and K+K'K^ decays based on a sample of 386 x lO^ BB 
pairs collected at the T(4S) resonance with the Belle detector at the KEKB energy-asymmetric 
e+e" collider. These decays are dominated by the b ^ s gluonic penguin transition and are sensi- 
tive to new CP-violating phases from physics beyond the standard model. One neutral B meson 
is fully reconstructed in one of the specified decay channels, and the flavor of the accompanying B 
meson is identified from its decay products. CP-violation parameters sinlcpl^ and Af for each of 
the decay modes are obtained from the asymmetries in the distributions of the proper-time inter- 
vals between the two B decays. We also perform an improved measurement of CP asymmetries 
in P° J/ipK^ decays using the same data sample. The same analysis procedure mentioned 
above yields sin 2(j)i = -1-0.652 =b 0.039(stat) ± 0.020(syst) , which serves as a reference point for the 
standard model, and Af = -FO.OlO ± 0.026(stat) ± 0.036(syst). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 



The flavor- changing b ^ s transition proceeds through loop penguin diagrams. Such 
loop diagrams play an important role in testing the standard model (SM) and new physics 
because particles beyond the SM can contribute via additional loop diagrams. CP violation 
in the h ^ s transition is especially sensitive to physics at a very high-energy scale [1]. 
Theoretical studies indicate that large deviations from the SM expectations are allowed for 
time-dependent CP asymmetries in meson decays [2]. Experimental investigations have 
recently been launched at the two B factories, each of which has produced more than 10^ BB 
pairs. The first measurement of the CP- violating asymmetry in decays [3], which 

are dominated by the b sss transition, by the Belle collaboration indicated deviation from 
the SM expectation [4]. Measurements with a larger data sample are required to confirm 
this difference. It is also essential to examine additional modes that are sensitive to the 
same b s penguin amplitude. In this spirit, experimental results based on a sample of 
275 X 10^ BB pairs using decay modes 5° 4)K^, t]'K% K^gK^sK^g, ir^vrO, f^K^ uK'^s, and 

K~ Kg [5] have already been reported [6, 7]. The combined result differs from the SM 
expectation by 2.4 standard deviations. Since measurements by the BaBar collaboration also 
yield a similar deviation [8, 9], the present world average differs from the SM expectation 
by 3.7 standard deviations. 

In the SM, CP violation arises from a single irreducible phase, the Kobayashi-Maskawa 
(KM) phase [10], in the weak-interaction quark-mixing matrix. In particular, the SM pre- 
dicts CP asymmetries in the time-dependent rates for B^ and B^ decays to a common CP 
eigenstate fcp [H]- In the decay chain T(45') B^B^ fcpftug, where one of the B 
mesons decays at time tc'p to a final state fcp and the other decays at time ttag to a final 
state /tag that distinguishes between 5° and 5°, the decay rate has a time dependence given 

by 

V{At) = — <^ 1 + q- SfSm{AmdAt) + AfCos{AmaAt) \. (1) 

Here Sf and Af are CP-violation parameters, tbo is the P° lifetime, Arud is the mass 
difference between the two P° mass eigenstates. At = tcp — ttagj and the 6-fiavor charge 
q = +1 (—1) when the tagging B meson is a P" (P*^)- To a good approximation, the 
SM predicts Sf = — ^/sin20i, where ^/ = 1) corresponds to CP-even (-odd) final 

states, and ^/ = for both b — > ccs and b — > sqq transitions. Therefore, a comparison of 
CP- violation parameters between b — > sqq and b — > ccs decays is an important test of the 
SM. 

Recent theoretical studies [12] find that P" — > (pK^, rj'K^ and KgKgKg have the smallest 
hadronic uncertainties among the modes hsted above. The effective sin20i values, sm2(j)f^, 
obtained from these decays arc expected to agree with sin 20i from the P° — > J/ ipK^ decay 
within 0.04. Larger deviations would indicate a new CP-violating phase beyond the SM. 
The other modes may be affected by a larger amount by the b ^ u transition that has a 
weak phase 03. Correspondingly, the SM predictions for the sin 20^^ values of these modes 
suffer larger uncertainties. 

Belle's previous measurements of CP violation in P° <t>K^, (f)K^, fi'Kg, KgKgK^, 
KgTi^, foKg, ujKg and K^K^Kg decays were based on a 253 fb^^ data sample containing 
275 X 10^ PP pairs. In this report, we describe improved measurements for these decays 
incorporating an additional 104 fb~^ data sample that contains 111 x 10^ PP pairs for a 
total of 386 X 10^ PP pairs. We also measure CP asymmetries for P° — > rj'K^ and rj'Kg 
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followed by Kg — > tt^tt^, which were not included in the previous analysis. 

Recent measurements of time-dependent CP asymmetries in decay modes governed by 
the b — >■ CCS transition by Belle [13, 14] and BaBar [15] have determined sin20i = +0.726 ± 
0.037 [9], where 5° ^ J/^K^;, J/i^K^, ip{2S)K^s, Xci^^s and r/ei^g decays are used. In 
this report, we describe improved measurements of CP-violation parameters Sf and Af in 
J/ipKg and J/ipK^ decays, which are the modes with the largest statistics and with 
the smallest theoretical uncertainties [16, 17], as a firm reference point for the SM. 

Among the 6 — > s modes listed above, all of the two-body final states are CP eigenstates 
with a CP eigenvalue ^/ = -1 r/'i^g, /T^tt^ and uK^) or ^/ = +1 {(f)Kl, rj'Kl and 

foK^). While the three body state K^K^K^s is a CP eigenstate with ^/ = +1 [18], the 
K~^K~Kg state is in general a mixture of both CP-even and -odd final states. Excluding 
K^K~ pairs that are consistent with a — > K^K^ decay from the P° — > K~^K~Kg 
sample, we find that the K'^K'Kg state is primarily CP-even; a measurement of the CP- 
even fraction using the isospin relation [19] with a 357 fb~^ data sample gives /+ = 
0.93±0.09(stat)±0.05(syst). The SM expectation for this mode is 5/ = -(2/+-1) sin20i. In 
this report, we define ^/ = 2/+ - 1 = +0.86±0.18(stat) ±0.09(syst) for the P° ^ K+R-R^s 
decay, and measure sin20f'^ = —^J^Sf. 

The decays P° — > (f)Rg and 0ir£ are combined in this analysis by redefining Sf as — ^/<S/ 
to take the opposite CP eigenvalues into account, and are collectively called (f)R^" . 

Likewise, CP asymmetries for rj'R^" or "P*^ J/ipR^" are obtained by combining 

the decays P° ^ r/'i^g and rj'Rl, or 5° ^ J/^Rs and J/^Rt 

At the KEKB energy-asymmetric e^e" (3.5 on 8.0 GeV) collider [20], the T(4S') is pro- 
duced with a Lorentz boost of (5^ — 0.425 nearly along the electron beamline {z). Since 
the P° and P° mesons are approximately at rest in the T(4S') center-of-mass system (cms). 
At can be determined from the displacement in z between the fcp and /tag decay vertices: 

~ {zcp - ^tag)/(/37c) = A^/(/?7c). 

The Belle detector is a large-sohd-angle magnetic spectrometer that consists of a siheon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect R^ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [21]. Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVD-I) were 
used for the first 140 fb~^ data sample (DS-I) that contains 152 x 10^ BB pairs, while a 

I. 5 cm radius beampipe, a 4-layer siheon detector (SVD-II) [22] and a small-cell inner drift 
chamber were used for the rest, a 217 fb~^ data sample (DS-II) that contains 234 x 10^ BB 
pairs. 

II. EVENT SELECTION, FLAVOR TAGGING AND VERTEX RECONSTRUC- 
TION 

A. Overview 

We reconstruct the following P° decay modes to measure CP asymmetries: B^ 4>Kg, 
(pRl, r]'R% r]'Rl, R^R^R^ R^g^^, f^R^g^ ooR^g and R+R-R% We exclude R+R- pairs 
that are consistent with a — > R'^R~ decay from the P° — > R'^R~R^g sample. The 
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intermediate meson states are reconstructed from the following decays: 7r° 77, Kg 
TT+TT", 7] 'j'j, p° — > TT+TT", uj 7r"'"7r~7r°, 1]' p^'~f OT r^TT+VT", /q TT+TT", and ^ K~ . 
In addition, i^^ 7r°7r° decays are used for — >• and r^'i^T^ decays, and rj — >• 7r+7r~7r° 
for the case 5° V^s (^s ^ Tr+Tr"). 

B. ct>K% and iiT+K-Kg 

Charged tracks reconstructed with the CDC for kaon and pion candidates, except for 
tracks from tt+tt^ decays, are required to originate from the interaction point (IP). 

We distinguish charged kaons from pions based on a kaon (pion) likelihood Ck{-k) derived 
from the TOF, ACC and dE/dx measurements in the CDC. 

Pairs of oppositely charged tracks that have an invariant mass within 0.015 GcV/c^ of 
the nominal mass are used to reconstruct tt+tt^ decays. The distance of closest 

approach of the candidate charged tracks to the IP in the plane perpendicular to the z axis 
is required to be larger than 0.02 cm for high momentum (> 1.5 GeV/c) Kg candidates 
and larger than 0.03 cm for those with momentum less than 1.5 GeV/c. The 7r"'"7r~ vertex 
is required to be displaced from the IP by a minimum transverse distance of 0.22 cm for 
high-momentum candidates and 0.08 cm for the remaining candidates. The mismatch in the 
z direction at the K^ vertex point for the tt+tt^ tracks must be less than 2.4 cm for high- 
momentum candidates and less than 1.8 cm for the remaining candidates. The direction of 
the pion pair momentum must also agree with the direction of the vertex point from the IP 
to within 0.03 rad for high-momentum candidates, and to within 0.1 rad for the remaining 
candidates. The resolution of the reconstructed Kg mass is 0.003 GeV/c^. 

Photons are identified as isolated ECL clusters that are not matched to any charged 
track. To select K^ tt^tt" decays, we reconstruct 7r° candidates from pairs of photons 
with > 0.05 GeV, where E^^ is the photon energy measured with the ECL. Photon pairs 
with an invariant mass between 0.08 and 0.15 GeV/c^ and a momentum above 0.1 GeV/c are 
used as 7r° candidates. Initially, the 7r° decay vertex is assumed to be the IP. An asymmetric 
mass window is used to take into account the lower tail of the mass distribution due to the 
distance between the IP and the true 7r° vertex. Candidate Kg 7r°7r° decays are required 
to have an invariant mass between 0.47 GeV/c^ and 0.52 GeV/c^, where we perform a fit 
with constraints on the Kg vertex and the tt^ masses to improve the tt^tt^ invariant mass 
resolution. We also require that the distance between the IP and the reconstructed Kg decay 
vertex be larger than —10 cm, where the positive direction is defined by the Kg momentum. 

Candidate K^K~ decays are required to have an invariant mass that is within 
0.01 GeV/c^ of the nominal meson mass. Since the (j) meson selection is effective 
in reducing background events, we impose only minimal kaon-identification requirements; 
TZk/tt = J^k/{J^k + J^tt) > 0.1 is required, where the kaon likelihood ratio TZk/tt has val- 
ues between (likely to be a pion) and 1 (likely to be a kaon). We use a more stringent 
kaon-identification requirement, TZk/w > 0.6, to select non-resonant K^K~ candidates for 
the decay — > K^K~Kg. We exclude K^K^ pairs with an invariant mass within 0.015 
GeV/c^ of the nominal meson mass to reduce the (p contribution to a negligible level. 
To remove Xco K~^K~, J/ip ^ K^K~ and — > K'^K~ decays, K^K^ pairs with 
an invariant mass within 0.015 GeV/c^ of the nominal masses of Xco and J/tJj or within 
0.01 GeV/c^ of the nominal mass are rejected. — > KgK~^ decays are also removed by 
rejecting KgK^ pairs with an invariant mass within 0.01 GeV/c^ of the nominal D'^ mass. 

For reconstructed B ^ fcp candidates, we identify B meson decays using the en- 
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ergy difference AE = — E^^^ and tfie beam-energy constrained mass Mbc = 

\/(-^bcam)^ ^ {pW^^y^ wliere i?bcam ^^c bcam energy in tlie cms, and E'j^^ and p^'' are 
tfie cms energy and momentum of tfie reconstructed B candidate, respectively. Tlie reso- 
lution of Mbc is about 0.003 GeV/c^. Because of the smallness of p^'^, the Mbc resolution 
is dominated by the beam-energy spread, which is common to all decay modes. The reso- 
lution in AE depends on the reconstructed decay mode. The AE resolution is 0.013 GeV 
for {K^s TT+TT-) and K+R-R^ The AE distribution for {R^ -k^-k^) has a 
tail toward lower AE due to 7 energy leakage in the ECL. The typical AE resolution for 
(PR's {Kg tt^tt") is 0.058 GeV for the main component and the typical width of the tail 
component is about 0.14 GeV. The B meson signal region is defined as \AE\ < 0.06 GeV 
for 5° ^ </)Kg (X^ ^ vr+TT-), -0.15 GeV < AE < 0.1 GeV for 5° ^fCg (K^ -> 7r°7r°), 
\AE\ < 0.04 GeV for 5° ^ R+R-R^, and 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ for all 
decays. 

The dominant background to the S° — > (f)Rg and R'^R'Rg decays comes from e+e~ — > 
v/u, dd, ss, or cc continuum events. Since these tend to be jet-like, while the signal events 
tend to be spherical, we use a set of variables that characterize the event topology to distin- 
guish between the two. We combine S±, 9t and modified Fox- Wolfram moments [23] into a 
Fisher discriminant where is the scalar sum of the transverse momenta of particles 
other than the reconstructed B candidate outside a 45° cone around the candidate (f) meson 
direction (the thrust axis of the B candidate for R^R^Rg decays) divided by the scalar sum 
of their total momenta, and 9t is the angle between the thrust axis of the B candidate and 
that of the other particles in the cms. We also use the angle of the reconstructed B candidate 
with respect to the beam direction in the cms (Ob)- We combine and cos 9b into a signal 
[background] likelihood variable, which is defined as £sig[bkg] = -Csigifakg] (-^) x ^sig[bkg](cos^B)- 
We impose requirements on the likelihood ratio TZs/h = ^sig/ (-^sig + ^hkg) to maximize the 
figure-of-merit (FoM) defined as N^^/ ^Jn^^P + iVbkg, where N^^^ {N^kg) represents the ex- 
pected number of signal (background) events in the signal region. We estimate using 
Monte Carlo (MC) events, while A^bkg is determined from events outside the signal region. 

We define two TZs/h regions for the decay S° — > (f)Rg {Rg — > tt+tt"). We require TZs/h > 
0.65 for the high-T^-g/b region. The requirement for the low-7^s/b region depends on the 
fiavor-tagging quality, r, which is described in Sec. II K. The threshold values range from 
0.1 (used for r > 0.875) to 0.35 (used for r < 0.25). For the S° ^ (pR^. {R^g 7r°7r°) 
candidates, the TZg/^, threshold values depend on r and range from 0.4 to 0.75, which are 
more stringent than those for the Rg tt+tt^ case. For the B^ R^R~Rg candidates, 
we require |cos^r| < 0.9 prior to the TZg/h requirement. The TZs/h threshold values range 
from 0.25 to 0.65. The TZs/h requirement reduces the continuum background by 65% for 
50 ^ (j)R"g {R°g TT+TT-), 92% for B^ R+R'R^g and 93% for 5° ^ (j)R^g {R°g ttV), 
retaining 91% of the signal for 5° ^ {R^ tt+tt"), 72% for 5° ^ R+R-R^ and 
78% for fiO ^ {Rl ^ ttOttO). 

We use events outside the signal region as well as a large MC sample to study the 
background components. The dominant background is from continuum. The contributions 
from BB events are small. We estimate the contamination of B^ — > R'^R''Rg and B^ — > 
/o^s (/o R+R-) decays in the B^ (pR^ sample from the Dalitz plot for B R+R'R 
candidates with a method that is described elsewhere [19]. The contamination of B^ — > 
R^R~Rg events in the B^ — > (pRg sample is 2.75±0.14%, which is taken into account in our 
signal yield extraction. The background fraction from the decay B° — > foRg (/o — > R'^R~), 
which has a CP eigenvalue opposite to (j)Rs, is found to be consistent with zero. The 
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FIG. 1: Distributions of (a) Mbc in the AE signal region, (b) AE in the Mbc signal region and (c) 
cos 6h in the Ai?-Mbc signal region for — > (pKg candidates. Solid curves show the fits to signal 
plus background distributions, and dashed curves show the background contributions. 

influence of the foKg background is treated as a source of systematic uncertainty. 

Figures 1(a), (b) and (c) show the distributions of Mbc in the AE signal region, AE 
in the Mbc signal region and cos6h in the Aii^-Mbc signal region for the reconstructed 
—>■ (pKg candidates. Here the helicity angle Oh is defined as the angle between the 
B meson momentum and the daughter momentum in the (p meson rest frame. The 
signal yield for the B^ — >■ (pK^ decay is determined from an unbinned three-dimensional 
maximum-likelihood fit to the AE-M\yc-cos6H distribution [24]. The fit region is defined as 
-0.12 GeV < AE < 0.25 GeV for the fsT^ ^ Tr+vr" channel, -0.25 GeV < AE < 0.25 GeV 
for the Kg —>■ tt^tt^ channel and Mbc > 5.2 GeV/c^ for both cases. The signal distribution for 
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FIG. 2: Distributions of (a) M^^ within the AE signal region, (b) AE within the Mbc signal region 
for K^K~ candidates. Solid curves show the fits to signal plus background distributions, 

and dashed curves show the background contributions. 

(Kg TT+TT^) is modeled with a Gaussian function (a sum of two Gaussian functions) 
for Mbc The (t)K^ {K^ -k^-k^) signal distribution is modeled with a smoothed 

histogram obtained from MC events. For the continuum background, we use the ARGUS 
parameterization [25] for Mbc and a linear function for AE. Finally, the cos Oh distribution 
for the signal (continuum) is modeled with a second-order polynomial and is 

determined from MC (events in the AE-M\^c sideband). The cosOh distribution for the 
no n- resonant K^K^Kg background is also determined from MC and is included in 

the fit, with a ratio between the non-resonant component and the 4>Kg signal fixed at the 
measured value. The fits yield a total of 180 ± 16(stat) B^ — > (pK^ events in the signal 
region. 

Figures 2(a) and (b) show distributions of Mbc in the AE signal region and AE in the 
Mbc signal region for the reconstructed B^ K^K~Kg candidates after flavor tagging and 
vertex reconstruction. The signal yield for the B^ — >■ K^K~Kg decay is determined from 
an unbinned two-dimensional maximum-likelihood fit to the Aii^-Mbc distribution in the fit 
region defined as -0.12 GeV < AE < 0.25 GeV and Mbc > 5.2 GeV/c^. The signal and 
background distributions are modeled in the same way as the B^ 0-^5 (-^s ~^ tt^tt^) 
case. The fit yields 536 ± 29(stat) B^ K~^K~Kg events in the signal region. 

C. S° cj)Kl 

Candidate K^K~ decays are selected with the criteria described above. We select 
candidates based on KLM and ECL information. There are three classes of candi- 
dates, which we refer to as KLM, ECL and KLM+ECL candidates. The KLM candidates 
are selected from hit clusters in the KLM that are not associated with either an ECL cluster 
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nor with a charged track. The requirements for the KLM candidates are the same as those 
used in the — > J/ipK^ selection for our previous sin20i measurement [13]. ECL candi- 
dates are selected from ECL clusters if there is no KLM candidate. We use a likelihood 
ratio [13], which is calculated from the following information: the distance between the ECL 
cluster and the closest extrapolated charged track position; the ECL cluster energy; Eg/E25, 
the ratio of energies summed in 3 x 3 and 5x5 arrays of CsI(Tl) crystals surrounding the 
crystal at the center of the shower; the ECL shower width and the invariant mass of the 
shower. The likelihood ratio is required to be greater than 0.69. A KLM+ECL candidate 
is an ECL cluster with cluster energy greater than 0.16 GeV that has an associated KLM 
cluster. Here we impose less stringent requirements than those for KLM candidates to select 
the cluster in the KLM detector. The likelihood ratio for the ECL cluster is required to 
be greater than 0.56. For all KLM, KLM+ECL and ECL candidates, we also require that 
the cosine of the angle between the direction and the direction of the missing momentum 
of the event in the laboratory frame be greater than 0.6. 

Since the energy of the is not measured, Afbc and AE cannot be calculated in the 
same way as for the other final states. Using the four-momentum of a reconstructed 4> 
candidate and the flight direction, we calculate the momentum of the K^^ candidate 
requiring AE — 0. We then calculate pg^^, the momentum of the B candidate in the cms, 
and define the B meson signal region as 0.2 GeV/c < p™^ < 0.5 GeV/c. We impose the 
requirement TZs/h > 0.80, which rejects 95.7% of the continuum background and 67.0% of 
backgrounds from B decays, while retaining 65.2% of signal events. Here TZg/h is based 
on the discriminating variables used for the S° — > (f)Kg decay and the number of tracks 
originating from the IP with a momentum above 0.1 GeV/c. We exclusively reconstruct 
and reject 5° ^ K+R-R^s (including ^ K+R- and /o ^ R+R-), <pR*^ {R*° R+n- 
or /sT^tt"), 07r°, 0r;, 5+ ^ and {R*+ R^s^+ or R+tt^) decays. If there is 

more than one candidate B^ — > (pRi^ decay in the signal region, priority is given to KLM 
candidates. If there still exist multiple candidates, we take the one with the R^ candidate 
closest to the expected R^ direction. 

We study the background components using a large MC sample as well as data taken 
with cms energy 60 MeV below the nominal T(4»S') mass (off-resonance data). The dominant 
background is from continuum. A MC study shows that background events from B decays 
are dominated by inclusive B — » (f)R^X decays that include B — > (f)R* decays. 

The signal yield is determined from an extended three-dimensional binned maximum- 
likelihood fit to the TZs/h-p's^^-r distribution in the fit region 0.8 < TZg/h < 1-0, GeV/c < 
^giis ^ Q g QqY/c and r > 0.25, where the total likehhood is a product of the likehhood for 
each of three variables. The S° — >■ (pR^ signal shape is obtained from MC events. Back- 
ground from BB pairs is also modeled with MC. We fix the ratio between the signal and 
the BB background based on known branching fractions and MC-determined reconstruction 
efficiencies with the R^ detection efficiency corrected from B^ J /ipR^ data. The uncer- 
tainty in the ratio is treated as a source of systematic error. The continuum background 
distribution is represented by a histogram obtained from MC events; we confirm that the 
function well describes both the off-resonance data and the events in a p^*^ sideband re- 
gion defined as 1.0 GeV/c < pg"'" < 1.6 GeV/c. The fit yields 78 ± 13 5° ^ (t)Rl events, 
where the error is statistical only. The result is in agreement with the expected B^ — > 0/^2 
signal yield (59 events) obtained from MC after applying the efficiency correction from the 
B^ J/ipR^ data. Figure 3(a) shows the 7^s/b distribution in the p'^^^^-r signal region. 
Figure 3(b) shows signal yields obtained for six p^*^ intervals separately. The yields agree 
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FIG. 3: (a) Distribution of TZ^/^ in the pg"^ signal region. The sohd histogram shows the fit to the 
signal plus background distribution, and the shaded histograms show the background contributions, 
(b) Background-subtracted p^*^ distribution for 'P^l candidates. The solid histogram shows 

the result of the three-dimensional fit to the 'R-s/h-p's^^-''" distribution. 

with the distribution obtained by the three-dimensional fit. 
D. B° ^ ry'Kg 

Candidate tt+tt^ and ti^ti^ decays are selected with the same criteria as those used 

for the — ^ 0-^5 decay. Charged pions from the 1], or r]' decay are selected from tracks 
originating from the IP. We reject kaon candidates by requiring TZk/tt < 0.9. Candidate 
photons from tt^ —>■ 77 decays are required to have > 0.05 GeV. The reconstructed 7r° 
candidate is required to satisfy 0.118 GeV/c^ < M^^ < 0.15 GeV/c^ and p^"?" > 0.1 GeV/c, 
where M^^ and p'^^ are the invariant mass and the momentum in the cms, respectively. 
Candidate photons from 77 — > 77 {rj' — > p°7) decays are required to have > 0.05 (0.1) GeV. 
The invariant mass of the photon pair is required to be between 0.5 and 0.57 GeV/c^ 
for the ?7 — > 77 decay. The vr+Tr^Tr" invariant mass is required to be between 0.535 and 
0.558 GeV/c^ for the rj Tr^Tr^Tr" decay, which is used only for the reconstruction of the 
j^o _^ rj'Kg {Kg —>■ TT+TT^) decay. A kinematic fit with an t] mass constraint is performed 
using the fitted vertex of the 7r"'"7r~ tracks from the rj' as the decay point. For rj' — > p°7 
decays, candidate p° mesons are reconstructed from pairs of vertex-constrained tt+tt^ tracks 
with invariant mass between 0.55 and 0.92 GeV/c^. The rj' —>■ rjn^n^ candidates are required 
to have a reconstructed mass between 0.94 and 0.97 GeV/c^ (0.95 and 0.966 GeV/c^) for 
the 77 — >• 77 (r] — >• n~^7i~n^) decay. Candidate rj' p°7 decays are required to have a 
reconstructed mass from 0.935 to 0.975 GeV/c^. 

The B meson signal region is defined as \AE\ < 0.06 GeV for B^ v'^s iv' ~^ 
p°7, K'^s vr+TT-), -0.1 GeV < AE < 0.08 GeV for S° ^ iv' ^ W^n', 7] 

77, K^s vr+TT-), -0.08 GeV < AE < 0.06 GeV for 5° ^ r^'Kg (r/' t/tt+tt-, t] 
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FIG. 4: Distributions of (a) Mbc within the AE signal region, (b) AE within the Mbc signal region 
for rj' candidates. Solid curves show the fits to signal plus background distributions, and 

dashed curves show the background contributions. 



tt+tt-ttO, Kl TT+TT-), -0.15 GeV < AE < 0.1 GeV for 5° ^ r^'iT^ (Kg n^n^), and 
5.27 GeV/ < Mbc < 5.29 GeV/ for all decays. The continuum suppression is based on the 
likelihood ratio TZg/h obtained from the same discriminating variables used for the 5° 4>Kg 
decay, except for the decay mode rj' pj {p ^ tt+tt") where cos9h is included. Here 9h is 
defined as the angle between the rj' meson momentum and the daughter vr"*" momentum in 
the p meson rest frame. The minimum T^-g/b requirement depends both on the decay mode 
and on the flavor-tagging quality, and ranges from (i.e., no requirement) to 0.4 for the de- 
cay B° 7]'K^s {K^s TT+TT-) and from 0.2 to 0.9 for the decay 5° r/'iTg (ifg ^ 7t°7t°). 
For the rj' — > p°7 mode, we also require |cos6't| < 0.9 prior to the T^g/b requirement. 
With these requirements, the continuum background in the — > Tj'K^ {K^ — > tt+tt") 
mode is reduced by 87% for rj' p°7, 58% for rj' rj-n^n^ {rj 77) and 31% for 
rj' rjn^n^ {rj tt+tt^tt^), while retaining 78% of the signal for rj' —>■ p°7, 94% for 
rj' rjTi^Ti^ {rj 77) and 97% for rj' —>■ //tt+tt" {rj — >■ tt+tt-tt^). The continuum background 
for the 5° rj'K^ {K^ tt°tt°) candidates is reduced by 90% (97%) while retaining 81% 
(54%) of signal events for rj' rj-n^n^ {pi)- 

We use events outside the signal region as well as a large MC sample to study the 
background components in v'^s- '^^^ dominant background is from continuum. In 

addition, according to MC simulation, there is a small (~ 3%) combinatorial background 
from BB events in B^ v'^s W ~^ P^l)- The contributions from BB events are smaller 
for other modes. The influence of these backgrounds is treated as a source of systematic 
uncertainty. 

Figure 4(a) shows the Mbc distribution for the reconstructed B^ v'^s candidates 
within the AE signal region after flavor tagging and vertex reconstruction, where all sub- 
decay modes are combined. The AE distribution for the B^ v'^s candidates within the 
Mbc signal region is shown in Fig. 4(b). The signal yields are determined from unbinned 
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two-dimensional maximum-likelihood fits to the Ai?-Afbc distributions. The fit region is de- 
fined as -0.25 GeV < AE < 0.25 GeV and Mbc > 5.2 GeV/c^. We perform the fit for each 
final state separately. The rj'Kg {Kg — 7i~^n~) signal distribution is modeled with a sum of 
two (three) Gaussian functions for Mbc i^E). The rj'Kg {Kg — > 7r°7r°) signal distribution is 
modeled with a smoothed histogram. For the continuum background, we use the ARGUS 
parameterization for Mbc and a linear function for AE. For the t]' fr/ mode, we include 
in the fits the BB background shape obtained from MC. The fits yield a total of 830 ± 35 
5° — > Tj'Kg events in the signal region, where the error is statistical only. 

E. S° rj'Kl 

Candidate rj' ?77r"'"7r~ (77 77) decays are selected with the same criteria as those 
used for the v'^s analysis. The K^ selection is adopted from the (l)Kl analysis, with 

a likelihood ratio optimized for the B^ — > rj'K^^ decay that is required to be greater than 
0.50 (0.40) for KLM+ECL (ECL) candidates. The best candidate is formed from the 77' 
candidate with the smallest value in its mass-constrained fit and the candidate whose 
measured direction is closest to the expected direction. The following exclusive modes 
are reconstructed and arc rejected: B t^'tt", ri'-n^, rj'ri, rj'Kg, Tj'K^. rj'K*^ {-^ K^tt^ 
or i^^TT^^), Tj'K*^, Tj' and rj'p^. The B meson signal region is defined as TZs/h > 0.8, 
0.2 GeV/c < p""^ < 0.5 GeV/c and r > 0.25 (0.5 for ECL candidates). 

The signal yield is determined from an extended three-dimensional maximum-likelihood 
fit to the TZs/h-p'^B^^-T distribution. The procedure to determine the signal and background 
distributions is the same as that for the B^ — > (pK^ decay. The fit yields 187 ± 18 5° — > 
•q'K^ events, where the error is statistical only. The result is in good agreement with the 
expected B^ — > r]'K^ signal yield (180 events) obtained from MC after applying the efficiency 
correction from the B^ J/ipK^ data. Figure 5(a) shows the TZs/h distribution in the p'^^^-r 
signal region. Figure 5(b) shows signal yields obtained for twelve p^** intervals separately. 
The yields agree with the distribution obtained by the three-dimensional fit. 

F. fiO ^ K%K%K% 

We reconstruct the 5° ^ K^K^K^ decay in the K^^' K-^' K-^' or K^'K^'Kf final 
state, where the tt+tt^ {-k^-k^) state from a K^ decay is denoted as K^^ i^s^)- Pahs 
of oppositely charged tracks with tt+tt" invariant mass within 0.012 GeV/c^ (~ 3a) of the 
nominal Kg mass are used to reconstruct Kg~ candidates. The tt'^tt" vertex is required to be 
displaced from the interaction point (IP) by a minimum transverse distance of 0.22 cm for Kg 
candidates with momentum greater than 1.5 GeV/c and 0.08 cm for those with momentum 
less than 1.5 GeV/c. The angle in the transverse plane between the Kg momentum vector 
and the direction defined by the Kg vertex and the IP should be less than 0.03 rad (0.1 rad) 
for the high (low) momentum candidates. The mismatch in the z direction at the Kg 
vertex point for the two charged pion tracks should be less than 2.4 cm (1.8 cm) for the 
high (low) momentum candidates. After two good Kg~ candidates have been found that 
satisfy the criteria given above, looser requirements are applied for the third Kg~ candidate. 
The requirement on the transverse direction matching is relaxed to 0.2 rad (0.4 rad for low 
momentum candidates), and the mismatch of the two charged pions in the z direction is 
required to be less than 5 cm (1 cm if both pions have hits in the SVD). We also require 
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FIG. 5: Distributions of (a) TZ^/h ™- the p™*^ signal region and (b) background-subtracted p'^^^ for 
Tj'K^ candidates. Solid histograms show the fits to signal plus background distributions, 
and the shaded histogram shows the background contributions. 



that the flight length in the plane perpendicular to the beam axis be less than 0.5 mm 
and the momentum be greater than 0.5 GeV/c. 

To select candidates, we reconstruct 7r° candidates from pairs of photons with > 
0.05 GeV. The reconstructed vr^ candidate is required to have an invariant mass between 
0.08 and 0.15 GeV/c^ and momentum above 0.1 GeV/c. Kg^ candidates are required to have 
an invariant mass between 0.47 and 0.52 GeV/c^, and a fit is performed with constraints 
on the Kg vertex and tt^ masses to improve the tt^tt" invariant mass resolution. The Kg^ 
candidate is combined with two good Kg~ candidates to reconstruct a meson. 

The meson signal region is defined as \^E\ < 0.10 GeV for B^ —>■ Kg~ Kg^ Kg~ , 
-0.15 GeV < AE < 0.10 GeV for 5° ^ K+~K+~K°g^, and 5.27 GeV/c^ < Mbc < 
5.29 GeV/c^ for both decays. To suppress the e~^e~ qq continuum background (g = 
u, d, s, c), we form the likelihood ratio TZs/h by combining likelihoods for two quanti- 
ties; a Fisher discriminant of modified Fox- Wolfram moments, and the cosine of the cms 
B^ flight direction. The requirement for TZg/h depends both on the decay mode and on the 
flavor-tagging quality; after applying all other cuts, this requirement rejects 94% of the qq 
background while retaining 75% of the signal. 

If both B^ Kg" Kg~ Kg~ and Kg~Kg~K^ candidates are found in the same event, we 
choose the B^ — > Kg^Kg^Kg~ candidate. If more than one B^ Kg~Kg~Kg~ candidate 
is found, we check for each of them the quality of the third Kg^ candidate, which is selected 
with looser requirements as described above. We choose the B^ — ^ Kg~Kg~Kg~ candidate 
in which the third Kg~ candidate satisfies the tight Kg~ selection requirements. If no B^ 
candidate is found with the tight requirements or more than one B^ candidate still remain, 
we select the one with the smallest value for ^{A.Mj^+-y' , where AM^+- is the difference 

between the reconstructed and nominal mass of Kg~. For multiple B^ Kg^Kg~K^ 
candidates, we select the Kg~Kg~ pair that has the smallest Y^{/S.Mj^+-Y value and the 
candidate with the minimum of the constrained fit. 
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FIG. 6: Distributions of (a) M^^ within the AE signal region, (b) AE within the Mbc signal region 
for — > KgKgKg candidates. Solid curves show the fits to signal plus background distributions, 
and dashed curves show the background contributions. 



We reject KgKgKg candidates if they are consistent with XcoKg {KgKg)Kg or 

— s> D^Kg — s> {KgKg)Kg decays, i.e. if one of the Kg pairs has an invariant mass within 
±2(7 of the Xco mass or mass, where a is the KgKg mass resolution. 

We use events outside the signal region as well as a large MC sample to study the 
background components. The dominant background is from continuum. The contamination 
of — i> XcoKg events in the B^ KgKgKg sample is small. The contributions from other 
BB events are negligibly small. The influence of these backgrounds is treated as a source of 
systematic uncertainty in the CP asymmetry measurement. Backgrounds from the decay 
B^ —>■ D^Kg are found to be negligible. 

Figure 6 shows the Mbc and AE distributions for the reconstructed B^ — > KgKgKg 
candidates after flavor tagging and vertex reconstruction. The signal yield is determined 
from an unbinned two-dimensional maximum-likelihood flt to the AE-M]^^ distribution. 
The Kg~ Kg~ Kg~ signal distribution is modeled with a Gaussian function (a sum of two 
Gaussian functions) for Mbc (A-E). For B^ —>■ Kg~Kg~K^ decay, the signal is modeled 
with a two-dimensional smoothed histogram obtained from MC events. For the continuum 
background, we use the ARGUS parameterization for Mbc and a linear function for AE. 
The fits after fiavor tagging and vertex reconstruction yield 88 ± 10 5° Kg~Kg~Kg~ 
events and 16 ± 6 S° ^ K^^-K^^-Kf events for a total of 105 ± 12 5° ^5^5^° events 
in the signal region, where the errors are statistical only. The obtained purity is 0.70 for 
the Kg~ Kg~ Kg" and 0.43 for the Kg" Kg~ Kg^ channels. Here the purity is defined as 
Ngig/Nev, where iVgig is the number of signal events in the signal region obtained by the fit, 
and A'^ev is the total number of events in the signal region. 
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FIG. 7: Distributions of (a) Mbc within the AE signal region, (b) AE within the Mbc signal region 
for — > K^TT^ candidates. Solid curves show the fits to signal plus background distributions, and 
dashed curves show the background contributions. 

G. ^ K^ttO 

Candidate Kg —>■ vr+vr" decays are selected with the same criteria as those used for the 
4>Kg decay, except that we use pairs of oppositely charged pions that have an invariant 
mass within 0.018 GeV/c^ of the nominal Kg mass. The 7r° selection criteria are the same 
as those used for the v'^s decay. 

The B meson signal region is defined as —0.15 GeV < AE < 0.1 GeV and 5.27 GeV/c^ < 
Mbc < 5.29 GeV/c^. The AE distribution for iT^vr^ has a tail toward lower AE. The AE 
resolution is 0.047 GeV for the main component. The width of the tail is about 0.1 GeV. The 
dominant background is from continuum. In addition, according to MC simulation, there 
is a small (~ 2%) contamination from other charmless rare B decays. We use extended 
modified Fox- Wolfram moments, which were applied for the selection of the B^ vr^vr^ 
decay [6], to form a Fisher discriminant JF. We then combine likelihoods for JF and cos 6b 
to obtain the event likelihood ratio TZs/h for continuum suppression. 

As described below, we include events that do not have B decay vertex information in 
our fit to obtain a better sensitivity for the CP-violation parameter Af. For events with 
and without vertex information, the high-7?.s/b region is defined as TZs/h > 0.8 and the low- 
TZs/h region as 0.45 < TZs/h < 0.8 for both DS-I and DS-II. After applying the high-7^s/b 
requirement, 95% of the continuum background is rejected and 62% of signal events remain. 
In the low-7?.s/b region, 84% of the continuum background is rejected and 24% of the signal 
remains. 

Figure 7(a) shows the Mbc distribution for the B^ —>■ K^n^ candidates within the AE 
signal region after flavor tagging and before vertex reconstruction. Also shown in Fig. 7(b) 
is the AE distribution within the Mbc signal region. The signal yield is determined from 
an unbinned two-dimensional maximum-likelihood fit to the AE-Mi^c distribution in the fit 
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region defined as 5.2 GeV/c^ < Mbc < 5.29 GeV/c^ and -0.3 GeV < AE < 0.3 GeV. 
The KgTi^ signal distribution is modeled with a smoothed histogram obtained from 

MC and calibrated with data using B~ — > D^tt~ {D^ — > K~tt^7t^). For the continuum 
background, we use the ARGUS parameterization for and a hnear function for AE. 
The B decay background distribution is represented by a smoothed histogram obtained 
from MC simulation. The fits yield 248 ±20 and 96 ±23 B° K'^n^ events in the high-7^s/b 
and low-T^-s/b signal regions, respectively, where the errors are statistical only. The same 
procedure after the vertex reconstruction yields a total of 106 ± 14 K^n^ events. 

H. SO^/oKg 

Candidate Kg -k^-k^ decays are selected with criteria that are slightly different from 
those used for the B^ decay so as to obtain the best sensitivity to CP violation in the 

B^ —>■ foKg decay. Pairs of oppositely charged tracks that have an invariant mass between 
0.484 GeV/c^ and 0.513 GeV/c^ are used to reconstruct Ks — > tt+tt" decays. The distance 
of closest approach of the candidate charged tracks to the IP in the plane perpendicular to 
z axis is required to be larger than 0.008 cm. The vr+vr" vertex is required to be displaced 
from the IP by a minimum transverse distance distance of 0.1 cm. The direction of the pion 
pair momentum must also agree with the direction of the vertex point from the IP to within 
0.03 rad. 

Pairs of oppositely charged pions that have invariant masses between 0.890 and 1.088 
GeV/c^ are used to reconstruct /o n^ii^ decays. Tracks that are identified as kaons 
{T^K/n > 0.7) or electrons are not used. We reject both K^n^ and KgTT~ combinations with 
an invariant mass within 0.02 GeV/c^ of the nominal charged D meson mass to remove 
background from — > Kgir^. 

The B meson signal region is defined as —0.03 GeV < AE < 0.06 GeV and 5.27 GeV/c^ < 
Mbc < 5.29 GeV/c^. The AE resolution is about 20 MeV. The dominant background is 
from continuum. The likelihood ratio TZg/h is obtained from cos^^, and cos9h, where 
the helicity angle 9h is defined as the angle between the B*^ meson momentum and the tt"*" 
momentum in the /q meson rest frame. The requirement for TZs/h depends on the flavor 
tagging r, and the threshold values range from 0.3 (used for r > 0.875) to 0.8 (used for 
r < 0.25). The continuum background is reduced by 93%, while retaining 72% of signal 
events with the requirement on TZg/i,- 

Figure 8(a) shows the Afbc distribution for the reconstructed S° — > foKg candidates 
within the AE signal region after flavor tagging and vertex reconstruction. The AE distri- 
bution for the B^ — > foKg candidates within the Mbc signal region is shown in Fig. 8(b). 
For the signal yield extraction, we first perform an unbinned two-dimensional maximum- 
likehhood fit to the A£'-Mbc distribution in the fit region defined as Mbc > 5.2 GeV/c^ 
and —0.12 GeV < AE < 0.3 GeV. The signal is modeled with a Gaussian function (a 
sum of two Gaussian functions) for Mbc (AE). For the continuum background, we use 
the ARGUS parameterization for Mbc and a linear function for AE. The fit yields the 
number of B^ ir'^ir^Kg events that have tt'^tt" invariant masses within the fo resonance 
region, which may include contributions from B^ — > p^Kg as well as non-resonant three- 
body B^ TT'^TT^Kg decays. To separate these peaking backgrounds from the B^ foKg 
decay, we perform another fit to the 7r~^7r~ invariant mass distribution for the events inside 
the A£'-Mbc signal region. We use Breit-Wigner functions for the B^ — > foKg signal, for the 
S° — > pKg background and for a possible resonance above the /o mass region, which is re- 
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FIG. 8: Distributions of (a) Mbc in the AE signal region, (b) AE in the Mbc signal region and (c) 
m^+^- in the A£^-Mbc signal region for fo^s candidates. Solid curves show the fits to signal 

plus background distributions. Dashed curves in (a) and (b) show the background contributions. 
In (c), each point shows a yield for 7r~^7r~Kg including foKg obtained from a fit to the 

A£^-Mbc distribution, the dotted line shows the B^ — > pKg, and the dashed line shows other quasi 
two-body decays as well as three-body B^ ir'^n^Kg decays. 



ferred to as /x(1300). The contributions from other resonant or non-resonant —>■ n~^7i~Kg 
decays are modeled with a threshold function. The combinatorial background is represented 
by the Mbc-AE" sideband and subtracted from the signal region distribution. The vr+vr" 
invariant mass distribution with the fit result is shown in Fig. 8(c). The fit yields 145 ± 16 
B° ^ foK^s events. 



I. B° 

Candidate Kg tt+tt^ decays are selected with criteria that are identical to those used 
for the 5° (pKg decay. Pions for the u vt+tt^tt^ decay are selected with the same 
criteria used for the r] n^n^n^ decay, except that we require p'j^^ > 0.35 GeV/c. The 
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7r"'"7r~7r° invariant mass M^^^ is required to be between 0.73 GeV/c^ and 0.84 GeV/c^. The 
B meson signal region is defined as —0.10 GeV < AE < 0.08 GeV and 5.27 GeV/c^ < 
Mbc < 5.29 GeV/c^. The AE resolution is 0.028 GeV. The dominant background is from 
continuum. The continuum suppression is based on the likelihood ratio TZg/h obtained from 
the same discriminating variables used for the — > (f)Kg decay plus the helicity angle 9h 
defined as the angle between the 5° meson momentum and the cross product of the tt"*" and 
7r~ momenta in the lu meson rest frame. We also require | cos^^l < 0.9 prior to the TZg/b 
requirement. We define two TZs/h regions. The TZg/h requirements depend on the flavor- 
tagging quality. The boundary between the high-7?.s/b regions and the low-T^-s/b regions is 
0.85 for all r values. The minimum TZs/h requirements range from 0.1 to 0.6 for the low-T^-g/b 
regions. The TZg/h and | cos^tI requirements reject 85% of the continuum background while 
retaining 84% of the signal. The contribution from BB events is negligibly small. 

Figures 9(a-c) show the Mbc distribution for the reconstructed B^ candidates 
within the AE signal region, the AE distribution within the Mbc signal region and the 
Mstt distribution within the AE-M\^c signal region, respectively, after flavor tagging and 
vertex reconstruction. The signal yield is determined from an unbinned three-dimensional 
maximum-likelihood flt to the AE-M\^c-^zv distribution in the flt region deflned as Mbc > 
5.2 GeV/c^, -0.12 GeV < AE < 0.25 GeV and 0.73 GeV/c^ < Mg^ < 0.84 GeV/c^. The 
B^ ujKg signal distribution is modeled with a sum of two (three) Gaussian functions for 
Mbc {AE and M^j^). For the continuum background, we use the ARGUS parameterization 
for Mbc, a linear function for AE and a second-order polynomial function plus three Gaussian 
functions for Ma^. The flt yields 68 ± 13 S° — > coKg events in the signal region. 

J. SO J/il^K% and J/il^Kl 

The reconstruction and selection criteria for B^ J/i/jK^ decays used in this mea- 
surement are the same as those in the previous publication, which are described in detail 
elsewhere [13]. We reconstruct J/ip candidates via their decays to {i = iJ,,e), and Kg 
candidates via Kg — > tt+tt" decays. The B meson signal region is deflned as |A£^| < 0.04 
GeV and 5.27 GeV/c^ < Mbc < 5.29 GeV/c^. 

Candidate J/iJj fi~^fi~ or e~^e~ decays for the B^ J/ipK^ mode are selected by 
requiring 3.05 GeV/c^ < M^^ < 3.13 GeV/c^ or 2.95 GeV/c^ < Mge < 3.13 GeV/c^, 
where M^^ (Mge) is the invariant mass of the (e^e^) pair. The momentum of the 

reconstructed J/V' candidate is required to be between 1.38 GeV/c and 2.00 GeV/c. The 
selection criteria for K^ candidates are identical to those in the B^ analysis, except 

that the K^ likelihood ratio for the ECL cluster is required to be greater than 0.25 for both 
KLM+ECL and ECL candidates. The B signal region is deflned as 0.2 GeV/c < p^'^ < 
0.45 GeV/c. 

Figure 10(a) shows the Mbc distribution for the reconstructed B^ J/i]jK^ candi- 
dates within the AE signal region after flavor tagging and vertex reconstruction. The AE 
distribution for the B^ — ^ J/tjjKg candidates within the Mbc signal region is shown in 
Fig. 10(b). The signal yield for the S° — > J/ipKg decay is determined from an unbinned 
two-dimensional maximum-likelihood flt to the A£^-Mbc distribution. The flt region is de- 
flned as \AE\ < 0.05 GeV and Mbc > 5.2 GeV/c^. The signal distribution is modeled with 
a Gaussian function (a sum of two Gaussian functions) for Mbc (AE). For the background, 
we use the ARGUS parameterization for Mbc and a linear function for AE. Figure 10(c) 
shows the distribution for the reconstructed S° — > J/ipK^^ candidates. The signal yield 
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FIG. 9: Distributions of (a) Mbc within the AE signal region, (b) AE within the My,c signal 
region and (c) m^+^-^o within the A£^-Mbc signal region for — > toKg candidates. Solid curves 
show the fits to signal plus background distributions, and dashed curves show the background 
contributions. 



for the — > J/ ipK^ decay is determined from a binned maximum-likelihood fit to the p^*^ 
distribution for each of KLM, KLM+ECL and ECL candidates separately. The fit region is 
defined as GeV/c < p^^^^ < 2 GeV/c. The shapes of the signal and background with J/ip 
are determined from the J / inclusive MC sample. Here background distributions with 
and without are treated separately to minimize the effect of an uncertainty in the 
detection efficiency in the MC simulation. The background shape for the case with a fake 
J/ ip meson is obtained from events in the sideband of the mass distribution. 
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FIG. 10: Distributions of (a) Mbc within the AE signal region and (b) AE within the Mbc signal 
region for J/tpKg candidates, and (c) p'j^^ for J/ipK'l candidates. 



K. Flavor Tagging 

The 6-flavor of the accompanying B meson is identified from inclusive properties of par- 
ticles that are not associated with the reconstructed 5° fcp decay. We use the same 
procedure as for our previous sin20i measurement [14]. The algorithm for flavor tagging is 
described in detail elsewhere [28]. We use two parameters, q and r, to represent the tagging 
information. The first, g, is defined in Eq. (1). The parameter r is an event-by-event, MC- 
determined flavor-tagging dilution factor that ranges from r = for no flavor discrimination 
to r = 1 for unambiguous flavor assignment. It is used only to sort data into six r intervals 
listed in Table I. The wrong tag fractions for the six r intervals, wi [l = 1,6), and differ- 
ences between 5° and decays, Awi, are determined from the data; we use the same values 
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TABLE I: Event fractions e;, wrong-tag fractions wi, wrong-tag fraction differences Awi, and 
average effective tagging efficiencies e^^^ = e;(l — 2wi)^ for eacfi r interval for DS-II. Errors for wi 
and Awi include both statistical and systematic uncertainties. The event fractions are obtained 



from J/ijjKg data. 


I 


r interval 


e« 




Awi 


^eff 


1 


0.000 - 0.250 


0.384 ±0.011 


0.467 ± 0.006 


±0.005 ± 0.007 


0.002 ± 0.001 


2 


0.250 - 0.500 


0.165 ±0.007 


0.324 ± 0.007 


-0.029 ± 0.009 


0.021 ± 0.002 


3 


0.500 - 0.625 


0.105 ±0.006 


0.223 ±0.010 


±0.019 ±0.011 


0.032 ± 0.003 


4 


0.625 - 0.750 


0.112 ±0.006 


0.160 ±0.011 


±0.008 ± 0.011 


0.052 ± 0.004 


5 


0.750 - 0.875 


0.089 ± 0.005 


0.101 ± 0.009 


-0.022 ± 0.010 


0.057 ± 0.004 


6 


0.875 - 1.000 


0.144 ±0.007 


0.020 ± 0.006 


±0.003 ± 0.006 


0.133 ± 0.007 



that were used for the sin20i measurement [14] for DS-I. Wrong tag fractions for DS-II are 
separately obtained with the same procedure and are listed in Table I. The total effective 
tagging efficiency for DS-II is determined to be e^Q = J2i=i ^/(l — 2w;)^ = 0.30 ± 0.01, where 
€i is the event fraction for each r interval determined from the J/ipKg data and is listed in 
Table I. The error includes both statistical and systematic uncertainties. We find that the 
wrong tag fractions for DS-II are slightly smaller than those for DS-I. As a result, the eggp 
value for DS-II is shghtly larger than that for DS-I (e^g = 0.287 ± 0.005). 

L. Vertex Reconstruction 

The vertex position for the fcp decay is reconstructed using charged tracks that have 
enough SVD hits: at least one layer with hits on both sides and at least one additional z 

hit in other layers for SVD-I, and at least two layers with hits on both sides for SVD-II. 
A constraint on the IP is also used with the selected tracks; the IP profile is convolved 
with the finite B flight length in the plane perpendicular to the z axis. The pions from 
decays are not used except in the analysis of — > -f^^TT^ and K^K^K^ decays. The typical 
vertex reconstruction efficiency and z resolution for B^ — >■ 4>K^ decays are 95% and 78 /im, 
respectively. Similar values are obtained for other fcp decays except for B^ ^s^^ 
K^K^K^ decays. 

The vertex for S° K'gH^ decays is reconstructed using the Kg trajectory and the IP 
constraint, where both pions from the decay are required to have enough SVD hits in 
the same way as for other fcp decays. The reconstruction efficiency depends both on the 
momentum and on the SVD geometry; the efficiency with SVD-II (32%) is significantly 
higher than that with SVD-I (23%) because of the larger outer radius and the additional 
layer. The typical z resolution of the vertex reconstructed with the is 93 //m for SVD-I 
and 110 yum for SVD-II. 

The vertex position for 5° KgKgKg decays is also obtained using Kg^ trajectories and 
a constraint on the IP. The reconstruction efficiency depends both on the Kg~ momentum 
and on the SVD geometry. The vertex efficiencies with SVD-I (SVD-II) are 79% (86%) 
for Kj-Kj'Kj' and 62% (74%) for K^'K^-Kf. The typical vertex resolution is about 
97 /xm (113 yum) for SVD-I (SVD-II) when two or three Kg~ candidates can be used. The 
resolution is worse when only one Kg~ can be used; the typical value is 152 /xm (168 /xm) 
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TABLE II: Estimated signal purities and signal yields A^gig in the signal region for each fcp mode 
that is used to measure CP asymmetries. The purity is defined as Ngig/Ney, where N^v is the total 
number of events in the signal region. Results for K^tt^ decays are obtained with samples 

after flavor tagging but before vertex reconstruction. Results for other decays are obtained after 
flavor tagging and vertex reconstruction. 



Mode 




purity 


N ■ 




-1 


0.57 


180 ± 16 


ct>Kl 


+1 


0.12 


78 ±13 




-1 


0.61 


830 ± 35 




+1 


0.26 


187 ± 18 




+1 


0.64 


105 ± 12 




-1 


0.25 


344 ± 30 


kK's 


+1 


0.47 


145 ± 16 




-1 


0.17 


68 ±13 


K+K-K^ 


+0.86 ±0.18 ±0.09 


0.55 


536 ± 29 




-1 


0.98 


5264 ± 73 




+1 


0.60 


4792 ± 105 



for SVD-I (SVD-II), which is comparable to the /tag vertex resolution. 

The /tag vertex determination with SVD-I remains unchanged from the previous pub- 
heation [6], and is described in detail elsewhere [29]; to minimize the effect of long-hved 
particles, secondary vertices from charmed hadrons and a small fraction of poorly recon- 
structed tracks, we adopt an iterative procedure in which the track that gives the largest 
contribution to the vertex is removed at each step until a good is obtained. The 
reconstruction efficiency was measured to be 93%. The typical z resolution is 140 /xm [13]. 

For SVD-II, we find that the same vertex reconstruction algorithm results in a larger 
outlier fraction when only one track remains after the iteration procedure. Therefore, in 
this case, we repeat the iteration procedure with a more stringent requirement on the SVD- 
II hit pattern; at least two of the three outer layers have hits on both sides. The resulting 
outlier fraction, which is described in Sec. Ill, is comparable to that for SVD-I, while the 
inefficiency caused by this change is small (2.5%). 

M. Summary of Signal Yields 

The signal yields that contribute to the determination of CP-violation parameters, A^^sig, 
for fcp decays are summarized in Table II. These yields are obtained after flavor 

tagging and vertex reconstruction for all modes except — > K^n^. As events with no 
vertex information reduce the statistical error on Ako-h-o significantly, we include them in 
the fit for the B^ K^tt^ decay. The signal purities are also listed in the table. The signal 
yields are all consistent with expected values that are obtained from previously measured 
branching fractions [9] and reconstruction efficiencies estimated from MC simulation studies. 
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III. RESULTS OF CP ASYMMETRY MEASUREMENTS 



We determine Sf and Af for each mode by performing an unbinned maximum- likelihood 
fit to the observed At distribution. The probability density function (PDF) expected for the 
signal distribution, Psig(^^; .4./, g, w;, Auj^), is given by Eq. (1) incorporating the effect 
of incorrect flavor assignment. The distribution is convolved with the proper-time interval 
resolution function i?sig(At), which takes into account the finite vertex resolution. 

For the decays 5° ^ cpKl, rj'K^s, v'^l, foK^s, ^^s, K+K'K'^, J/il^K^g and J/tl^K^, 
we use fiavor-specific B decays governed by semileptonic or hadronic b ^ c transitions 
to determine the resolution function. We perform a simultaneous multiparameter fit to 
these high-statistics control samples to obtain the resolution function parameters, wrong- 
tag fractions (Section II K), Arrid, rB+ and r^o. We use the same resolution function used for 
the sin 20i measurement for DS-I [14]. For DS-II, the following modifications are introduced: 
a sum of two Gaussian functions is used to model the resolution of the fcp vertex while a 
single Gaussian function is used for DS-I; a sum of two Gaussian functions is used to model 
the resolution of the tag-side vertex obtained with one track and the IP constraint, while 
a single Gaussian function is used for DS-I. These modifications are needed to account for 
differences between SVD-I and SVD-II, as well as different background conditions in DS-I 
and DS-II. We test the resolution parameterization using MC events on which we overlay 
beam-related background taken from data. A fit to the MC sample yields correct values for 
all parameters. 

For the S° K^n^ decay, we use the resolution function described above with addi- 
tional parameters that rescale vertex errors. The rcscaling function depends on the detector 
configuration (SVD-I or SVD-II), SVD hit patterns of charged pions from the Kg decay, 
and Kg decay vertex position in the plane perpendicular to the beam axis. The parameters 
in the rescaling function are determined from a fit to the At distribution of B° — > J/ijjKg 
data. Here only the Kg and the IP constraint are used for the vertex reconstruction, the 
lifetime is fixed at the world average value, and 6-flavor tagging information is not used so 
that the expected PDF is an exponential function convolved with the resolution function. 

We check the resulting resolution function by also reconstructing the vertex with leptons 
from J/t/j decays and the IP constraint. We find that the distribution of the distance 
between the vertex positions obtained with the two methods is well represented by the 
resolution function convolved with the well-known resolution for the J/ip vertex. Finally, 
we also perform a fit to the J/ipKg sample with 6-flavor information and obtain 

Sj/^ko — +0.73±0.08(stat) and Aj/^ko = +0.01±0.04(stat), which are in good agreement 
with our measurement using leptons from J/ip decays, which will be described later. A 
separate fit to the same sample with tbo as a free parameter yields tbo = 1.55 ± 0.05(stat) 
ps, which is consistent with the world average value. Thus, we conclude that the vertex 
resolution for the — > KgH^ decay is well understood. 

For B^ — > KgKgKg candidates, we use the same resolution function that is used for the 
B^ — > KgTi^ decay if only one Kg~ is available for the vertex reconstruction. For events 
with n (= 2 or 3) Kg~ trajectories used in the vertexing, we adopt a function defined as 
[^(crj,)]^/" to rescale the vertex error cr^. Here V{az) is the aforementioned rescaling function 
for the case that only one Kg~ is available. We find from MC simulation that the resolution 
is well described by this form for the rescaling function. 
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We determine the following likelihood for each event: 



Pi = (I-/01) y[/signig(At')i?sig(At,- AO 

+ (1 - /sig)nkg(Ai')i?bkg(Aii - At') d{At') 

+ /oiPoi(At,), (2) 

where Poi(At) is a broad Gaussian function that represents an outlier component with a 
small fraction f^i [14]. The width of the outlier component for DS-I is determined to be 
(39li3) Ps; the fractions of the outher components are (2.lto;8) ^ 10~^ for events with the 
/tag vertex reconstructed with more than one track, and (3.llo;i) x 10"^ for the case only one 
track is used. Here the errors include both statistical and systematic errors. Corresponding 
values for DS-11 are (35l?i) ps, (3.6^?;?) x 10"^ and {l-StoD x lO^^. The signal probability 
/sig depends on the r region and is calculated on an event-by-event basis as a function of 
for the 5° J/^Kl decay, and 7^,/b for the S° (pK^ and r]'Kl decays, AE, 
Mbc and cos 6h for the B^* (p'Kg decay, AE, A'/^c and M^^r for the wKg decays, and AE 
and Mbc for the other modes. A PDF for background events, 'Pbkg(At), is modeled as a sum 
of exponential and prompt components, and is convolved with a sum of two Gaussians i^bkg- 
Parameters in 7^bkg(At) and i?bkg for continuum background are determined by a fit to the 
At distribution for events outside the A£'-Mbc signal region except for the S° — > (f)K^ and 
r]'Kl decays. For the B'^ 4>Kl and rj'Kl decays, we use sideband events to obtain 
the parameters. Parameters in 'Pbkg(At) and i?bkg for BB background events in B^ v'^s^ 
KgTT^, (pKl and r/'i^^ decays are determined from MC simulation. 

We fix tbo and Am^ at their world average values [30]. We assume no CP asymmetry in 
the background At distributions and possible CP asymmetries in the B decay backgrounds 
are treated as sources of systematic error. In order to reduce the statistical error on we 
include events without vertex information in the analysis of B^ — > Kg-K^. The likelihood in 
this case is obtained by integrating Eq. (2) over Atj. 

The only free parameters in the final fits are Sf and Af, which are determined by max- 
imizing the likelihood function L = Pi{Ati; Sf, Af) where the product is over all events. 
Table 111 summarizes the fit results of Sf and Af. We define the raw asymmetry in each 
At bin by {Nq=^i — iVg=_i)/(iVg=_)_i -|- iVq=_i), where A^g=_|_i(_i) is the number of observed 
candidates with q — 1). Figures 11-14 show the raw asymmetries for each decay mode 
in two regions of the fiavor-tagging parameter r [31]. Figures 15-17 also show At distribu- 
tions and asymmetries for B^ — > (pK'^, rj'K^ and J/ipK^ decays after subtracting background 
contributions, where the sign of each At measurement for the final states with is inverted 
to combine final states with and K^. 

Tables IV and V list the systematic errors on <S/ and Af, respectively. The total system- 
atic errors are obtained by adding each contribution in quadrature, and are smaller than the 
statistical errors for all b ^ s modes. 

To determine the systematic error that arises from uncertainties in the vertex recon- 
struction, the track and vertex selection criteria are varied to search for possible systematic 
biases. Small biases in the Az measurement are observed in e'^e~ — > jJi^ ^~ and other con- 
trol samples. Systematic errors are estimated by applying special correction functions to 
account for the observed biases, repeating the fit, and comparing the obtained values with 
the nominal results. The systematic error due to the IP constraint in the vertex reconstruc- 
tion is estimated by varying (±10//m) the smearing used to account for the B flight length. 
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TABLE III: Results of the fits to the At distributions. The first error is statistical and the second 
error is systematic. 



Mode 


SM expectation for 5/ 


Sf 


^/ 




+sin20i 


+0.44 ± 0.27 ± 0.05 


+0.14 + 0.17 + 0.07 


4>Kl 


+sin 201 


+0.19 ±0.32 


+0.12 + 0.20 




—sin 2(j)i 


-1.54 + 0.59 


+0.38 + 0.36 




+sin 201 


+0.62 + 0.12 + 0.04 


-0.04 + 0.08 + 0.06 




+sin20i 


+0.60 ± 0.14 


-0.04 + 0.09 


rfKl 


—sin 201 


-0.73 ± 0.29 


-0.02 + 0.18 




—sin 201 


-0.58 + 0.36 + 0.08 


+0.50 + 0.23 + 0.06 




+sin 201 


+0.22 + 0.47 + 0.08 


+0.11 + 0.18 + 0.08 




—sin 201 


-0.47 + 0.36 + 0.08 


-0.23 + 0.23 + 0.13 




+sin 201 


+0.95 + 0.53l{^;}? 


+0.19 + 0.39 + 0.13 




— (2 /_ — Ijsiii 2c;i 


-0.52 ± O.IG ± 0.03 


-0.06 ± 0.11 ± 0.07 




+sin 201 


+0.652 + 0.039 + 0.020 


+0.010 + 0.026 + 0.036 




+sin 201 


+0.668 + 0.047 


-0.021 + 0.034 




—sin 201 


-0.619 + 0.069 


+0.049 + 0.039 



Systematic errors due to imperfect SVD alignment are determined from MC samples that 
have artificial misalignment efi'ects to reproduce impact-parameter resolutions observed in 
data. 

Systematic errors due to uncertainties in the wrong tag fractions are studied by varying 

the wrong tag fraction individually for each r region. Systematic errors due to uncertainties 
in the resolution function are also estimated by varying each resolution parameter obtained 
from data (MC) by ±lcr (±2cr), repeating the fit and adding each variation in quadrature. 
Each physics parameter such as r^o and Am^ is also varied by its error. A possible fit bias 
is examined by fitting a large number of MC events. 

Systematic errors from uncertainties in the background fractions and in the background 
At shape are estimated by varying each background parameter obtained from data (MC) 
by ±1(7 (±2a). 

The PDF's for S° — > 0ir° and "q'Kl assume no correlation among T^-g/b, P^b^ and r. To 
estimate systematic errors due to possible correlations between TZ^/h and p^'', we repeat 
a fit to obtain CP parameters using signal fractions determined by the 7^s/b distribution 
for each p^*^ region separately. The difference from our nominal result is included in the 
systematic error. Systematic errors due to other possible correlations are estimated from 
events in the sideband, events with r < 0.25 and off-resonance data. 

Additional sources of systematic errors are considered for B decay backgrounds that are 
neglected in the PDF. We consider uncertainties both in their fractions and CP asymmetries. 
The effect of backgrounds from K^R-K^ and foK^g (/o ^ K^K') in the S° (pK^g 
sample is considered. Uncertainties from B — > (f)K* and other rare B decay backgrounds in 
the B^ 4>Kl sample are also taken into account. For the B^ v'^s sample, non-resonant 
B decay backgrounds are studied using events in the sideband of the reconstructed rj' mass 
distribution. Effects of possible CP asymmetries in B decay backgrounds for K^tt'^ and 
foKg are evaluated. The peaking background fraction in the S° — > foKg sample depends 
on the functions used to fit to the tt'^tt" invariant mass distribution. The systematic errors 
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FIG. 11: Raw asymmetry in each At bin with < r < 0.5 (top) and with 0.5 < r < 1.0 (bottom) 
for (a) B° -> (pK'^, (b) B° (pKl, (c) S° ^ rj'K^ and (d) S° ^ r/'K£ decays. The sohd 
curves show the results of the unbinned maximum-hkehhood fits. The dashed curves show the 
SM expectation with our measurement of CP-violation parameters for the J/ipK^ mode 

(sin 201 = +0.652 and Aj = +0.010). 



due to the uncertainties of the masses and widths of the resonances used in the fit are also 
included. The width of /o as well as the mass and the width of /x(1300) are varied by their 
errors. The effect of possible interference between resonant and non-resonant amplitudes, 
which is neglected in the nominal analysis, is also evaluated. We perform a fit to the tt+tt" 
distribution of a MC sample generated with interfering amplitudes and phases for B — >■ Kmr 
decays measured from data [19]. The observed difference in the signal yield from the true 
value is taken into account in the systematic error determination. We also repeat the fit to 
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FIG. 12: Raw asymmetry in each At bin with < r < 0.5 (top) and with 0.5 < r < 1.0 
(bottom) for (a) 5° ^ K^s^s^^s (b) ^° ^ K^s'^^, (c) /o^g, (d) 5° ^ toK^g and (e) 



— > K^K decays. The sohd curves show the results of the unbinned maximum-hkehhood 
fits. The dashed curves show the SM expectation with our measurement of CP-violation parameters 
for the B^ J/ipK^ mode (sin2(/)i = +0.652 and Af = +0.010). 



the At distribution ignoring the contribution of the peaking background. The differences in 
Sf and Af from our nominal results are included in the systematic error. 

Finally, we investigate the effects of interference between CKM-favored and CKM- 
suppressed B ^ D transitions in the /tag final state [32]. A small correction to the PDF for 
the signal distribution arises from the interference. We estimate the size of the correction 
using the D*^l^v sample. We then generate MC pseudoexperiments and make an 

ensemble test to obtain systematic biases in Sf and Af. In general, we find effects on Sf 
are negligibly small, while there are sizable possible shifts in Af. 
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FIG. 13: At distributions in (a) — > J/t/jK'^, (b) J/iJjK^, and raw asymmetries in (c) 

B^ J /ijjK^ and (d) B^ — > J/jpK^ with < r < 0.5. The curves show the results of the unbinned 
maximum-likehhood fits. 



TABLE IV: Summary of the systematic errors on Sf. 

















K+K-K^s 




Vertex reconstruction 


0.02 


0.02 


0.04 


0.01 


0.02 


0.02 


0.02 


0.015 


Flavor tagging 


0.01 


< 0.01 


0.01 


0.01 


0.01 


0.04 


< 0.01 


0.006 


Resolution function 


0.02 


0.02 


0.05 


0.05 


0.02 


0.05 


0.02 


0.005 


Physics parameters 


< 0.01 


< 0.01 


< 0.01 


0.01 


0.01 


0.01 


< 0.01 


0.001 


Possible fit bias 


0.01 


0.02 


0.03 


0.01 


0.02 


+0.01 
-0.09 


0.01 


0.007 


Background fraction 


0.01 


0.02 


0.02 


0.01 


0.02 


0.09 


0.02 


0.005 


Background At shape 


0.03 


< 0.01 


0.02 


0.06 


0.07 


0.01 


0.01 


0.006 


Tag-side interference 


< 0.01 


< 0.01 


0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


0.003 


Total 


0.05 


0.04 


0.08 


0.08 


0.08 


+0.12 
-0.15 


0.03 


0.020 



TABLE V: Summary of the systematic errors ov^ Af. 









K%K%K% 




hK% 




K+K-K^s 




Vertex reconstruction 


0.04 


0.04 


0.04 


0.04 


0.05 


0.07 


0.04 


0.030 


Flavor tagging 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.007 


Resolution function 


0.01 


< 0.01 


0.02 


< 0.01 


0.02 


0.03 


0.01 


0.001 


Physics parameters 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


0.001 


Possible fit bias 


0.01 


0.01 


0.02 


< 0.01 


0.01 


0.02 


0.01 


0.005 


Background fraction 


0.01 


0.01 


0.03 


0.01 


0.01 


0.10 


0.01 


0.004 


Background At shape 


0.02 


< 0.01 


0.01 


0.03 


0.09 


0.01 


< 0.01 


0.005 


Tag-side interference 


0.05 


0.04 


0.05 


0.06 


0.07 


0.04 


0.06 


0.017 


Total 


0.07 


0.06 


0.06 


0.08 


0.13 


0.13 


0.07 


0.036 
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FIG. 14: At distributions in (a) — > J/i/jK^, (b) i?^ ^ J/ipK^, and raw asymmetries in (c) 
—>■ J ftp Kg and (d) B^ — > J/ipK^ with 0.5 < r < 1. The curves show the results of the unbinned 
maximum-likehhood fits. 
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FIG. 15: Background-subtracted At distribution (top) and asymmetry in each At bin (bottom) 
with 0.5 < r < 1.0 for (j)K^ . The result of the unbinned maximum-hkehhood fit is also 

shown. The dashed curve in the bottom figure shows the SM expectation with our measurement 
of CP-violation parameters for the — s- J/ipK^ mode {sm2cl)i = -1-0.652 and Aj = -1-0.010). 
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FIG. 16: Background-subtracted At distributions and asymmetry in each At bin with 0.5 < r < 1.0 
for Tj'K^. The result of the unbinned maximum-hkehhood fit is also shown. The dashed curve 

in the bottom figure shows the SM expectation with our measurement of CP-violation parameters 
for the S° J/i/jK'^ mode (sin 201 = -^0.652 and Af = -^0.010). 
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FIG. 17: Background-subtracted At distributions and asymmetry in each At bin with 0.5 < r < 1.0 
for — > J/ijjK^ . The result of the unbinned maximum-hkehhood fit is also shown. 
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TABLE VI: Results of the sin 205 measurements. The first errors are statistical, the second errors 
are systematic, and the third error for the K'^K~Kg mode arises from the uncertainty in the 

CP-cvvii fracdou. 



Mode sin 2(f)f 

(j)K^ +0.44 ±0.27 ±0.05 

ri'K° +0.62 + 0.12 + 0.04 

K^K^K^ +0.58 + 0.36 ± 0.08 

K^TT^ +0.22 + 0.47 ± 0.08 

foK^ +0.47 + 0.36 ± 0.08 

wETg +0.95 + 0.53to:ll 

K+K-Kl +0.60 + 0.18 ±0.04l[i:^ 

J/i^K^ +0.652 + 0.039 + 0.020 



Various crosschecks of the measurements are performed. We reconstruct charged B meson 
decays that are the counterparts of the ^ fcp decays and apply the same fit procedure. 
All results for the Sf term are consistent with no CP asymmetry, as expected. Lifetime 
measurements are also performed for the fcp modes and the corresponding charged B 
decay modes. The fits yield r^o and Tb+ values consistent with the world average values. 
MC pseudocxpcriments are generated for each decay mode to perform ensemble tests. We 
find that the statistical errors obtained in our measurements are all consistent with the 
expectations from the ensemble tests. 

The results in this report are consistent with those in our previous publications [6, 7, 14] 
within statistical fluctuations and supersede them. Among our new results, the largest 
difference from the previous measurement is observed in the B^ — > KgKgKg decay. A fit to 
the 253 fb"^ data sample, which contains the entire DS-I and a part of DS-II and were used in 
the previous publication, yields Sf — +1.05±0.64(stat) and Af — +0.51±0.30(stat), where 
a small change from the previous measurement is due to an improvement in the selection 
of Kg candidates. A fit to an additional 104 fb~^ data sample alone yields Sf — —2.95 ± 
0.53(stat) and Aj = +0.64 ± 0.36(stat). Prom MC pseudoexperiments, the probability that 
the significance of difference is larger than the observed difference is estimated to be 1.0%. To 
check if this arises due to a difference in the SVDl and SVD2 detectors, we perform separate 
fits to DS-1 and DS-11. We obtain Sf = -0.38 ± 0.82(stat) and Af = +0.16 ± 0.43(stat) 
for DS-I and Sf = -0.92 ± 0.49(stat) and Af = +0.72 ± 0.28(stat) for DS-II, which are 
consistent to each other. As all the other checks mentioned above also yield results consistent 
with expectations, we conclude that the observed change in the CP-violation parameters 
for the B^ K^K^Kg mode is due to a statistical fiuctuation. 

Table VI summarizes the sin 20^^ determination based on our Sf measurements. For each 
mode, the first error shown in the table is statistical and the second error is systematic. For 
the P° K+K-K^ decay, the SM prediction is given by Sf = -(2/+ - l)sin20f . The 
third error is an additional systematic error arising from the uncertainty of the CP-even 
fraction. The results for each individual decay mode are consistent with sin20i obtained 
from the P° — > J/ipK^ decay within one standard deviation. 



34 



IV. SUMMARY 



We have performed improved measurements of CP-violation parameters sin202^ and Af 
for 5° ^ (l>K^, rj'K^, K^s^s^^s, ^s^r^, foK'^, uK^s and K+R-R^ decays. These charmless 
decays are dominated by 6 — > s flavor-changing neutral currents and are sensitive to possible 
new CP- violating phases. 

We have also measured CP asymmetries in — > J/ipR^ decays using the same data 
sample. The same analysis procedure as that used for the b s modes yields sin20i = 
+0.652 ±0.039(stat)±0.020(syst), which serves as a SM reference point, and Af = +0.010± 
0.026(stat) ± 0.036(syst). 

We do not see any significant deviation between the results for each b ^ s mode and 
those for P° — > J/ipR^. Since some models of new physics predict such effects, our results 
can be used to constrain these models. However, many models predict smaller deviations 
which we cannot rule out with the current experimental uncertainty. Therefore, further 
measurements with larger data samples are required in order to search for new, beyond the 
SM, CP-violating phases in the b ^ s transition. 

Acknowledgments 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the National Institute of Informatics for valuable computing and Super-SINET network 
support. We acknowledge support from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan and the Japan Society for the Promotion of Science; the Australian 
Research Council and the Australian Department of Education, Science and Training; the 
National Science Foundation of China under contract No. 10175071; the Department of 
Science and Technology of India; the BK21 program of the Ministry of Education of Korea 
and the CHEP SRC program of the Korea Science and Engineering Foundation; the Polish 
State Committee for Scientific Research under contract No. 2P03B 01324; the Ministry of 
Science and Technology of the Russian Federation; the Ministry of Higher Education, Science 
and Technology of the Republic of Slovenia; the Swiss National Science Foundation; the 
National Science Council and the Ministry of Education of Taiwan; and the U.S. Department 
of Energy. 



[1] A. G. Akeroyd et ai, "Physics at super B factory," hep-ex/0406071 and references therein. 
[2] See for example, 

Y. Grossman and M. P. Worah, Phys. Lett. B 395, 241 (1997); D. London and A. Soni, Phys. 
Lett. B 407, 61 (1997); T. Moroi, Phys. Lett. B 493, 366 (2000); D. Chang, A. Masicro and 
H. Murayama, Phys. Rev. D 67, 075013 (2003); S. Baek, T. Goto, Y. Okada and K. Okumura, 
Phys. Rev. D 64, 095001 (2001). 
[3] Throughout this paper, the inclusion of the charge conjugate decay mode is implied unless 
otherwise stated. 

[4] BcUc Collaboration, K. Abe et ai, Phys. Rev. Lett. 91, 261602 (2003). 

[5] Throughout this paper, 0, /o and u refer to 0(1020), /o(980) and a;(782), respectively. 



35 



[6] Belle Collaboration, K. F. Chen et al, hep-ex/0504023. 
[7] Belle Collaboration, K. Sumisawa et al, hep-cx/0503023. 

[8] BaBar Collaboration, B. Aubert et al., Phys. Rev. Lett. 94, 191802 (2005); liep-ex/0503011; 

hep-ex/0502019; hep-ex/0502013; hep-ex/0408095. 
[9] Heavy Flavor Averaging Group, J. Alexander et al, hep-ex/0412073. 
[10] M. Kobayashi and T. Maskawa, Prog. Thcor. Phys. 49, 652 (1973). 

[11] A. B. Carter and A. I. Sanda, Phys. Rev. D 23, 1567 (1981); I. I. Bigi and A. I. Sanda, Nucl. 
Phys. B193, 85 (1981). 

[12] M. Beneke and M. Neubert, Nucl. Phys. B675, 333 (2003); M. Beneke, hep-ph/0505075; 

S. Mishima, talk given at the second joint workshop on a super S-factory, April 2005, 

Hawaii, http://www.phys.liawaii.edu/supcrb/2005/slides/physics4/mishima/; H.-Y. Cheng, 

C.-K. Chua and A. Soni, hcp-ph/0502235; hep-ph/0506268. 
[13] Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 87, 091802 (2001); Phys. Rev. D 66, 

032007 (2002). 

[14] Belle Collaboration, K. Abe et al, Phys. Rev. D 71, 072003 (2005). 

[15] BaBar Collaboration. B. Aubert et al, Phys. Rev. Lett. 89, 201802 (2002); hep-ex/0408127. 

[16] H. Boos, T. Manuel and J. Renter, Phys. Rev. D 70, 036006 (2004). 

[17] D. Atwood and G. Hiller, hep-ph/0307251. 

[18] T. Gershon and M. Hazumi, Phys. Lett. B 596, 163 (2004). 

[19] Belle Collaboration, A. Garmash et al, Phys. Rev. D 69, 012001 (2004); Belle Collaboration, 

K. Abe et al, Bellc-CONF-0410 (2004). 
[20] S. Kurokawa and E. Kikutani, Nucl. lustrum. Methods Phys. Res., Sect. A 499, 1 (2003), and 

other papers included in this volume. 
[21] Belle Collaboration, A. Abashian et al, Nucl. lustrum. Methods Phys. Res., Sect. A 479, 117 

(2002). 

[22] Y. Ushiroda (Belle SVD2 Group), Nucl. lustrum. Methods Phys. Res., Sect. A 511, 6 (2003). 

[23] Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 87, 101801 (2001). 

[24] The signal yield for the 4>K% (-^,9 ^ tt^tt^) decay is determined from an unbinncd 

two-dimensional maximum-likelihood fit to the Ai?-Mb(, distribution because of the limited 
statistics. We include cos Oh in the Tls/h calculation for this mode to obtain the best sensitivity. 
The <pKg {Kg — > 7r°7r°) candidates are included in Fig. 1(a) and (b), but are not in (c). 

[25] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 241, 278 (1990). 

[26] Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 91, 261801 (2003). 

[27] T. Skwarnicki, Ph.D. Thesis, Institute for Nuclear Physics, Krakow 1986; DESY Internal 
Report, DESY F31-86-02 (1986). 

[28] H. Kakuno et al, Nucl. lustrum. Methods Phys. Res., Sect. A 533, 516 (2004). 

[29] H. Tajima et al, Nucl. lustrum. Methods Phys. Res., Sect. A 533, 370 (2004). 

[30] Particle Data Group, S. Eidclman et al. Particle Listings in the 2005 Review of Particle 
Physics, http: / / www.slac.stanford.edu/xorg/hfag/ osc/PDG_2005_draft6/. 

[31] While the numbers of signal events in the two r regions are similar, the effective tagging 
efficiency is much larger and the background dilution is smaller in the region 0.5 < r < 
1.0. Note that these projections onto the At axis do not take into account event-by-event 
information (such as the signal fraction, the wrong tag fraction and the vertex resolution), 
which is used in the unbinned maximum- likelihood fit. 

[32] O. Long, M. Baak, R. N. Cahn and D. Kirkby, Phys. Rev. D 68, 034010 (2003). 

[33] G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 (1998). 



36 



